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ABSTRACT 

We study Modified Chaplygin Gas (MCG) as a candidate for dark energy and 
predict the values of parameters of the gas for a physically viable cosmological model. 
The equation of state of MCG (p = Bp — -4r) involves three parameters: B : A and 
a. The permitted values of these parameters are determined with the help of dimen- 
sionless age parameter (H t ) and H(z) — z Data. Specifically we study the allowed 
ranges of values of B parameter in terms of a and A s (A s is defined in terms of the 
constants in the theory) . We explore the constraints of the parameters in Cold Dark 
Matter(CDM) model and UDME (Unified Dark Matter Energy) model respectively. 
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1 INTRODUCTION 

Recent cosmological observations, such as high redshift sur- 
veys of SNe la (Perlmutter et al. 1997a; Permutter et al. 
1997b; Riess et al. 1998; Tonry et al. 2003), CMBR (Mel- 
chiorri et al. 2000; Lange 2001; Jaffe et al. 2001; Netterfield 
et al. 2002; Halverson et al. 2002), WMAP (Briddle et al. 
2003; Bennet et al. 2003; Hinshaw et al. 2003; Kogut et al. 
2003; Spergel et al. 2003) predict that our present universe is 
passing through an accelerated phase of expansion preceeded 
by a period of deceleration. It is known that the ordinary 
matter and fields of the standard model are not sufficient to 
accommodate the present phase of acceleration (preceede by 
deceleration) . Consequently a modification of matter sector 
of the Einstein Gravity is essential to incorporate the recent 
prediction from observational cosmology. The notion of new 
type of matter has come up which must have negative pres- 
sure. Recent astronomical data when interpreted in the con- 
text of Big Bang Model have provided some interesting infor- 
mation about the composition of the universe. The analysis 
reveals that our universe is spatially flat and consists of 70 
percent dark energy with negative pressure, remaining 30 
percent dust matter (cold dark matter plus baryons), and 
negligible radiation. It has been predicted that the dark en- 
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ergy may be responsible for the present acceleration of our 
universe. 

The most simple candidate for these uniformly dis- 
tributed (i.e unclustered) dark energy is considered to be in 
the form of vacuum energy density or cosmological constant 
(A). The model with cosmological constant is entangled with 
(i) fine tuning problem (present amount of the dark energy 
is so small compared with fundamental scale) and (ii) co- 
incidence problem (dark energy density is comparable with 
critical density today). Alternatively the other choices are 
(i) a light homogeneous scalar field <f>, whose effective po- 
tential V(4>) leads to an accelerated phase at a later stage of 
the universe (Caldwell et al. 1998; Sahani et al. 2000), (ii) a 
X-matter component, which is characterized by an equation 
of state p = uip, where, — 1 uj < (Peebles and Ratra 
2002), (iii) effects from extra dimensions (Sahni 2002; Lui 
2002), (iv) an exotic fluid, Chaplygin Gas etc.. 

The motivation of the paper is to obtain cosmological 
model constraining Chaplygin Gas taking into account ob- 
servational facts. Chaplygin Gas was first introduced in aero- 
dynamics in 1904. Recently it has been shown that Chaply- 
gin Gas may be useful for describing dark energy because of 
its negative pressure. Although it has positive energy den- 
sity it carries a negative pressure for which it is referred to 
as exotic fluid. In the context of string theory the Chaply- 
gin gas emerges from the dynamics of a generalized d-brane 
in a (d+1,1) spacetime. It can be described by a complex 
scalar field which is obtained from a generalized Born-Infeld 
action. The equation of state is given by 
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p 



A 
P 



(1) 



where A is a positive constant, p and p are pressure and 
density respectively. Subsequently a modified form of the 
equation of state (Bento et al. 2002; Bilic, Tupper 2001) of 
the form 

P=-— (2) 

p a 

with < a ^ 1 which is known as Generalized Chaplygin 
Gas (GCG). It has two free parameter A (positive), a. In 
GCG model, at low energy density, the fluid pressure is neg- 
ative and constant while at high energy density it behaves 
almost like pressure less fluid. Thus it smoothly interpolates 
between a non-relativistic matter phase in the past and a 
negative pressure dark energy regime at late times. 

Recently a modified form of GCG has been considered 
in cosmology (Liu and Li 2005). The modified Chaplygin gas 
(MCG) is more general and contains three free parameters. 
The idea is to interpolate states of standard fluids at high 
pressures and at high energy densities to a constant nega- 
tive pressure at low energy densities (Debnath et al. 2004). 
In addition it covers whole aspects of GCG. This model ac- 
commodates consistent (i) Gravitational lensing test (Silva 
et al. 2003; Dev et al. 2004), (ii) Gamma-ray bursts (Berto- 
lami and Silva 2006) . The equation of state for this Modified 
Chaplygin Gas is given by 



p = Bp- 



A_ 

P a 



(3) 



where A, B, a are arbitrary constants with ^ a ^ 1. 
As there are three free parameters, unlike GCG, we look 
for a suitable range of B parameter for MCG for a viable 
cosmological model accommodating the observational evi- 
dences. The parameters are determined by (i) Considering 
a dimensionless age parameter H t (Dev et al. 2002) and 
(ii) H(z) - z Data analysis (Wu and Yu 2006). 
We perform the following analysis: 

Case 1. The age parameter (H t ) is dimensionless and 
a constant irrespective of the model we are considering. For 
simplicity we chose its standard value to be 0.95 (ignoring 
error). Imposing this constant age parameter we determine 
the effective ranges of values of the free parameters in this 
model. As the parameters have some preferred range of val- 
ues, one can ultimately constrain one of them, in particular 
the matter part B. 

Case 2. Using (H(z) — z) Data we further verify the 
validity of the constraints on the parameters obtained in 
case 1. We use Hubble parameter vs redshift relation given 
in Table 1. The \ 2 minimization technique has been used in 
this process. There are 9 data points of H{z) at redshift z 
used to constrain the MCG model. 

We investigate both CDM and UDME (Unified Dark 
Matter Energy) models in next sections. UDME model refers 
to the model in which the Modified Chaplygin gas(MCG) 
represent dark matter and dark energy as a whole, where 
the total energy density comprises of radiation, baryon and 
MCG energy density. In case of Cold dark matter (CDM) 
model the constituents of our universe are radiation, CDM 
and MCG. 

This paper is organized as follow : In section II we 
present the relevant field equations and introduce Hubble 



Table 1. 



z Data 


H(z) 


a 


0.09 


69 


± 12.0 


0.17 


83 


± 8.3 


0.27 


70 


± 14.0 


0.40 


87 


± 17.4 


0.88 


117 


± 23.4 


1.30 


168 


± 13.4 


1.43 


177 


± 14.2 


1.53 


140 


± 14.0 


1.75 


202 


± 40.4 



parameter and deceleration parameter respectively. In sec- 
tion III and IV we explore values of B parameters using the 
age parameter and the observed H(z) — z data respectively. 
In section V we checked viability of the models from union 
compilation data. In section VI we summerise the result. 



2 FIELD EQUATIONS, DECELERATION 
PARAMETER 

Let us now consider the Fricdmann-Robertson- Walker 
(FRW)line element (c=l) 

„2 



ds 2 = —dt 2 + a 2 (t) 



1 - kr 2 



+ r 2 (d6 2 + sin 2 6 dcj> 2 ) 



(4) 



where k = 0, ±1 is the curvature parameter in the spatial 
section a(t) is the scale factor of the universe, r, 6, <f> are the 
comoving co-ordinates. The energy conservation equation is 
given by 



dp 
~dl 



+ 3H(p + p) = 0, 



(5) 



where p, p, H are pressure, energy density, Hubble param- 
eter respectively. Using eq. (3) in eq. (5) we obtain the ex- 
pression for the energy density of MCG with the scale factor 
of the universe, which is given by, 



P = 



A _l 

l + B + a 



(6) 



where C is an arbitrary constant and we denote (1 + B)(l + 
a) = n. Equation (6) can be rewritten as 



P= Po 



A s + 



As 



(7) 



i 

l+z 



where z is redshift parameter, As = j^g +i i ^ 
and we chose ao = 1 for convenience. It reduces to GCG 
model when we set B — in the above equation. The Fried- 
mann's equation becomes 



H(z) = H [Q. r0 {l + z) 4 + 0,0(1 + zf + 
(1 - O,o - n j0 )[(A s + (1 - A s )(l + z ) 3n )TT^]]i 
The above equation can be rewritten in terms of a as 

H{a) - H [— + — 



(8) 



+(1 - Oo - Ojo)(A s + ■ 



A s 



V 



(9) 
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Figure 1. CDM, Dotted line for a = 0.0, Dashed line for a = 
0.2,Thin line for a = 0.39 




Figure 2. UDME,Thin line for a = O.O.Dotted line for a = 
0.5, Dashed line for a = 1.0 



where j = m for CDM model and j = b for UDME model. 
The above reduces to GCG model when one sets B = 0. The 
deceleration parameter (go = — (ft )t ) at the present time 
can be written as 



90 



Oo + f Oo + (1 + B)(l - Q j0 - Oo)(l - A s ) 



QjO + ^Cg + O 



-1(10) 



where j = m for CDM model and j = b for UDME 
model. The deceleration parameter can be estimated both 
in CDM and UDME model. For a flat universe we have 
fljo + Ocgo + Oo = 1 which will be used to measure the pa- 
rameters in the next section. In the above ficg represents 
the present day Modified Chaplygin gas energy density, Oo 
is the present energy density of either Cold dark matter (in 
CDM model) or baryonic energy density (in UDME model) 
and Oo represents the present radiation energy density of 
our universe. 



3 AGE OF OUR UNIVERSE AS A 
CONSTRAINING TOOL 

Using the definition of the age parameter 



to = 





da 




aH(a) 



(11) 



where 



y^-j and H(a) is given by eqn. (9). The predicted 
age of the universe in MCG model becomes 



tin 

with 



da 



af (a, Q joMrO, A s , B, a) 



f(a,tt j() ,Q r0 ,A s ,B,a) = 



H(a) 
H • 



(12) 



(13) 



We consider here from experimental facts the value H t = 
0.95 . Although it has some error limit in both sides, we take 
this value as standard. From the constancy of this parame- 
ter, we derive constraints on the parameters of the theory. 
For a given value of alpha we plot variation of A s with B. 
We note the following : 

Fig.l : shows variation of B with A s for a — 0, 0.20, 0.39 
by dotted, dashed and thin lines respectively in CDM model. 
It is evident that as the value of As approaches 1 (0.97 to 
1) for ^ a ^ 0.39 the B parameter picks up positive 
values with a maximum 0.20. 

Fig. 2: shows variation of B with A 3 for different a 
a=0,0.5 and 1 with thin, dotted and dashed lines respec- 
tively in UDME. In this case as the value of As is increased 



from 0.7 to 1 it is evident that the B parameter picks up 
positive value up to a maximum 1.02 for ^ a ^ 1 in 
UDME model. 

The set of curves shown in fig.l and fig. 2 is useful to de- 
termine the range of values of B for both CDM and UDME 
models respectively. We note that in CDM model B lies 
between to 0.20, where as in UDME model B lies be- 
tween to 1.02. Moreover, in CDM model B is positive 
when < a < 0.39 and A s between 0.97 to 1. In UDME we 
note that B is positive for < a < 1 and A s between 0.7 to 
1. 



4 H(Z) Z DATA AS CONSTRAINING TOOL 

For a flat universe containing only radiation, cold matter 
(or baryon) and the MCG, the Friedmann equation can be 
expressed as 

H 2 {Ho,A s ,B,a,z) = H 2 f(A s ,B,a,z) 



where, 
f(A a ,B,a,z) 



[0-0(1 + zf + £1,0 (1 + zf + 



(1 - Oo - Oo)(A + (1 - A a )(l + z) 



,l+al2 



(14) 



(15) 



with j = m for CDM model and j = b for UDME model. 
The best fit values for model parameters A s , B, a and Ho 
can be determined by minimizing as follows 

\ - (H(H ,A s ,B,a,z)~H obs (z)) 2 



X 2 (H ,A s ,B,a,z) 



(16) 



Since we are interested in determining the model pa- 
rameters, Ho is not an important parameter here. So we 
marginalize over Ho to evaluate the probability distribution 
function for A 3 , B, a as 



L(A s ,B,a) 



dH P(Ho)e X p( 



-X 2 (H ,A s ,B,a,z) 



(17) 



where P(Ho) is the prior distribution function for the 
present Hubble constant. We consider Gaussian priors Ho = 
72 ±8. For — la level calculation the limit of integration will 
be from 64 to 72 and for la level calculation the limit of 
integration becomes 72 to 80. Minimizing \ 2 determines the 
maximum L(A S , B , a) value. We determine the maximum 
value of the function L(A s ,B,a) by plotting the function 
with any of its parameter keeping other two fixed. As a result 
we get the maximum values of the parameters A s , B and a. 
Consequently a relation between B and A s for various a can 
be established. 
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Figure 3. CDM, Dashed line for 3a level, Dotted line for 2a 
level, Thin line for la level, for a = 0.01 




Figure 6. UDME, Dotted pair of lines represent boundary of 
l<T,Thin line pair represent 2a, Dashed pair for 3a level, for a = 




-0 . 6 



Figure 4. CDM, Thin line for 3a level, Dotted line for 2a 
level, Dashed line for la level, for a = 0.5 



B 




Figure 7. UDME, Dashed pair of represent boundary of lcr,Thin 
line pair represent 2a, Dotted pair for 3a level, for a = 0.5 



In CDM model, variation of B with A 3 for a = 
0.01, 0.5, 0.99 at -la, -2a and -3a level respectively are 
shown in figs. (3)-(5). We note that as the value of A$ tends 
to 1 we see that the B parameter picks up positive values (i) 
up to 1.07 (fig.3), (ii) upto 0.62 (fig.4), (iii) upto 0.36 (fig.5) 
in accordance with the H(z) — z data. Thus as a increases, 
B decreases. In UDME model variation of B with A 3 for 
a = 0, 0.5, 1 at Tier, ±2<r and ±3cr level respectively are 
shown in figs. (6)-(8). We note that as the value of As tends 
to 1 we see that the B parameter picks up positive values (i) 
upto 1.35 (fig.6), (ii) upto 0.84 (fig.7), (iii) upto 0.58 (fig.8) 
in accordance with the H(z) — z data. Thus as a increases, B 
decreases but compared to CDM model B parameter values 
in UDME is more for a given a and A 3 . 

Thus in CDM the range for B is to 1.07 and in Unified 
Dark Energy model the range is to 1.35 upto 99.7 percent 
confidence level. 

B 

/ 

0.2 /' 




In CDM B is positive only when A s is within 0.76 to 
1 for a to 1 and in UDME B is positive (so, permissible) 
only when A s is within 0.57 to 1 for a lying between to 1. 



5 VIABILITY OF THE MODEL FROM UNION 
COMPILATION OF SUPERNOVA 
MAGNITUDES AND REDSHIFT DATA 

We have already found the best-fit values of the parameters 
of MCG in our models from H(z)-z data part. For CDM 
model the best-fit values of the parameters are A s = .99, 
B = .01, a = .01 and for UDME model it is A s = .8,5 = .06, 
a — .11. In order to check its validity we used these best- 
fit values to find supernovae magnitudes at different red- 
shift for the two models of our consideration. From there we 
have drawn supernovae magnitudes vs redshift curve using 
those best-fit values of the models. We compared these su- 
pernovae magnitudes vs. redshift curve of those two models 



B 




Figure 5. CDM, Thin line for 3a level, Dashed line for 2a 
level, Dotted line for 3a level for a = 0.99 



Figure 8. UDME, Dashed pair of represent boundary of lcr,Thin 
line pair represent 2a, Dotted pair for 3a level, for a = 1.0 
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0.25 0.5 0.75 1 1.25 1.5 

Figure 9. fJ.(z) vs. Z curves for CDM model and union compila- 
tion data (in the figure m stands for fi(z)) 




0.25 0.5 0.75 



1.25 1.5 



Figure 10. fi(z) vs. Z curves for UDME model and union com- 
pilation data (in the figure m stands for fJ,(z)). 



to the original curve of union compilation data (Kowalaski 
et al.) (between those two parameters). Fig. 9 shows a plot 
of u(z) vs. z obtained from the CDM model (the continu- 
ous line) along with that obtained from union compilation 
data (the dots). The same curve is drawn for the UDME 
model in Fig. 10 (continuous line for UDME model and dots 
for union compilation data). As one can see from these two 
curves, both the CDM and UDME models are in excellent 
agreement with union compilation data. 



6 DISCUSSION 

Cosmological models with Modified Chaplygin Gas pre- 
sented here, contain three different parameters, A s , B, a. 
From the age constancy we determine permissible range of 
values of B parameter. In sec. 3 we plot B vs. A s for differ- 
ent values of a in figs. (1) and (2). The figures are plotted 
for both positive and negative values of B. In the case of 
CDM we note that B can pick up positive values upto 0.2 
for .97 < A s < 1, < a < 0.39. However, for UDME model 
we note that B can pick up positive values upto 1.02 for 
0.7 < A 3 < 1, < a ^ 1. 

In sec. (4), using the Hubble parameter vs. redshift data 
we obtain the constraints on B using chi-square minimiza- 
tion technique. To obtain a viable cosmology with MCG we 
restrict to positive values of B. The constraints on B are : 
(i) s: B sC 1.07 for 0.76 < A s < 1, and < a < 1 in CDM 



and (ii) < B < 1.35 for 0.56 < A s < 1, < a 1 for 
UDME. For UDME model the range of values of B is found 
to be more than that of CDM. If the age constant parameter 
is decreased then we note that the values of B permitted by 
CDM and UDME model are in agreement with that found 
by chi-square minimization of the observed H(z) vs z data. 
Consequently the limiting value of the age of our universe is 
pushed to lower values (t < 13.6 Billion years). 

The best-fit values of the parameters obtained here 
for CDM and UDME models are in agreement with union 
compilation data. We note that the best-fit values of our 
models are A s = 0.99,5 = 0.01, a = 0.01 for CDM and 
A s = 0.8,5 = 0.06,a = 0.11 for UDME model. 



ACKNOWLEDGMENTS 

PT and SG would like to thank IUCAA Reference Centre, 
Physics Department, N.B.U for extending the facilities of 
research work. SG would like to thank University of North 
Bengal for awarding Junior Research Fellowship. 



REFERENCES 

Riess A.G. et al.,1998,Astron.J., 116, 1009 

Perlmutter S. et al,1997,Bull. Am. Astron. Soc, 29,1351 

Perlmutter S. et al, 1997, Astr ophysics. J . 517,565 

Tonry J.L. et al,2003 |astro-ph/0305008"1 

Melchiorri A. et al.,2000,Astrophysics J.Lett.,536,L63 

Lange A.E. et al.,2001,Phys Rev. D, 63,042001 

Jaffe A.H. et al.,2001,Phys Rev. Lett. ,86,3475 

Netterfield C.B. et al.,2002,Astrophysics,J., 571,604 

Halverson N.W. et al., 2002, Astrophysics. J. ,568,38 

Briddle S. et al.,2003 Science, 299, 1532 

Bennet C. et al. ,2003|astro-ph703l)2207 ] 

Hinshaw G. et al. ,20 03|astro-ph/0302217l 

Kogut A. et al.,2003 |astro-ph/0302213| 

Spergel D.N. et al.,2003. ast"ro-ph/0302209l 

Caldwell R.R. et al.,1998,Phys.Rev.Lett.80,1582 

Sahani T.D. et al., 2000, Phys.Rev.Lett.85, 1162 

Sahani T.D. et al.,2000,Phys.Rev.Lett.85,1162 



Peebles P.J.E. and Ratra B.,2002,astro-ph/0207347 
Jain D. et al.,2001,astro-ph/010555 



Sahni V., Shtanov Y.,2002,astro-ph/0202346 0208823, 
Lue Arthur,2002 |hep-th/0208169| 

Kamenshchik A.,Moschella U. and Pasquier V.,2001,Phys. 

Lett B 511 265 
Bento M.C. et al.,2002,Phys.Dev. Dev. D 66 043507 
Bilic N. et al.,2001,Phys. Lett. B 535 17 
Liu Dao-Jun and Li Xin- 

zhou,2005 |arXiv:astro-ph/0501115| 

Debnath Ujjal. et al.,2004 |arXiv:gr-qc/041l6"l5"] 

Silva P.T et al.,2003,Astron. Astrophysics. 599 829 

Dev A. et al.,2004,Astron.Astrophys. 417 847 

Bertolami O. et al.,2006,Mon. Not. R. Astron, Soc. 365 

1149 

Dev Abha et al.,2002 p[st"ro-ph/0209379l 
Wu Puxun et al.,2006 ]arXiv:gr-qc/0612055| 
M.Kowalaski et al.,2008,arXiv:astro-ph/0804.4142 



